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INTRODUCTION _ 

The  question  of  shock  front  rise  time  is  fundamental  to  the  physics  of 
dynamic  high  pressure  and  high  strain  rate  effects  in  condensed  media.  The 
rise  time  question  is  currently  THE  basic  question  in  shock  wave  physics. 
Applications  range  from  a  better  understanding  of  the  shock  wave  initiation 
of  explosives  to  the  design  of  radically  new  sensors  for  impact  fuzing. 


In  a  classic  paper  Cowan  and  Homig  presented  a  study  of  the  optical 
reflectivity  of  a  shock  front  in  a  gas.  The  theory  presented  was  for  the 
linear  case  of  small  reflectivity  and  constant  angle  of  incidence  (within 
the  shock  front).  The  question  of  different  shock  structures  (i.e. ,  non¬ 
constant  density  gradients)  was  also  considered.  In  a  subsequent  paper^ 
plook  and  Hornig  attempted  to  measure  the  reflectivity  of  shock  fronts  in 
!a  number  of  transparent  liquids  including  water.  For  the  water  experiments 
the  reflectivity  was  too  small  to  be  measured. 

^  This  paper  reports  on  successful  shock  front  optical  reflectivity 
Jtheory  and  experiment  for  water  at  5.8  kbar  (0.58  GPa),  and  complimentary 
{shock  induced  electrical  polarization  experiment  and  theory  for  the 
'pressure  range  20  kbar  to  100  kbar.  The  LCWSL  Universite  de  Poitiers 
Iprogram  has  experiments  performed  in  France  (funded  by  the  Commission 
Rationale  Recherche  Sclent if ique),  while  the  ILIR  program  within  the  LCVTSL 
{funds  the  theoretical  work. 

I  ,  ^ 

!  ">The  consistent  interpretation  of  the  reflectivity  and  shock  polariza¬ 
tion  results  is  (a)  a  shock  front  rise  time  of  approximately  10~12  sec 
throughout  the  pressure  range  of  0  to  100  kbar,  (b)  a  transition  in  the 
{source  of  the  dominant  polarization  voltage  to  the  alignment  of  unbonded 
{dipole  moments  at  approximately  35  kbar,  and  (c)  in  the  unbonded  region 
‘the  average  contributing  dipole  is  perfectly  aligned  in  the  direction  of 
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The  shock  front  |e|lectivity  theory  and  experiment  have  already  been 


reported  on  elsewhere 
here. 


an/1 


and  consequently  that  work  will  only  be  summarized 


jreflectivity  of  a  (laser)  optical  pulse  from  a  propagating  5.8  khar  shock 


wave  front  In  water.  The  experiments  represent  the  first  time  that 
optically  polarized  beams  (both  parallel  and  perpendicular  to  the  plane  of 
incidence)  were  employed,  and  the  theoryr.  again-  for  the  first  time, 
[utilized  an  Index  of  refract  ion,  and-ang  3  e  of  incidence  which  varied  within 
the  thickness  of  the  shock  front. 


The  theory  can  be  summarized  via 


Reflectivity 
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where  the  integration  in  Eq.  (lc)  is  to  a  depth  z  within  the  shock  front 
thickness*  E  is  the  electric  field  intensity,  n  and  0  are  local  index  of 
tefractlon  and  angle  of  incidence  respectively,  Ya  an<*  Yg  are  fractional 
^optical  phase  angles,  and  the  subscript  ”o,v  denotes  preshock  values  while 
the  subscript  "r?>  denotes  a  reflected  parameter.  The  -  and  +  in  Eq.  (lb) 
denote  parallel  and  perpendicular  to  the  plane  of  incidence  respectively. 


Both  constant  gradient^ (i.e.  -g—  «  const)  and  hyperbolic  tangent 

jmodels  for  the  index  of  refraction  were  used  in  Eqs.  (1)  with  essentially 
the  same  analytic  result.  Aside  from  being  the  first  accurate  observation 
’and  analysis  of  shock  front  optical  reflectivity  for  condensed  media  to 
(appear  in  the  literature,  the  results  shown  in  Figs.  I  and  2  are  Important 
jin  that  they  provide  a  shock  front  thickness  upper  bound  of  approximately 


{400  A  at  5.8  kbar. 
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.  The  400  A  upper  bound  Is  extremely  Important  in ""that It  pins  the  shock) 
polarization  results;  to  now  we  have  only  been  successful  in  obtaining  valid 
shock  polarization  data  down  to  20  kbar,  j 

SHOCK  POLARIZATION 

The  shock  polarization  results,  including  theoretical  predictions  at 
[low  pressures  *  are  shown  in  Pig.  3.  A  modified  "rigid  four-point-charge” 
‘model®  has  been  used  in  analyzing  the  data,  and  in  making  the  low  pressure 
predictions.  _ 


In  homogeneous  media  computer  molecular  dynamics  (CMD)  calculations^ 
leach  four-point  ^charge  molecule  interacts  vita  neighboring  molecules  via 
Jan  electrostatic  potential  (as  well  as  via  an  oxygen~oxygen  interatomic 
jpotential) .  In  the  derivation  outlined  below  a  single  molecule  is  con¬ 
sidered  in  the  presence  of  a  shock  front  induced  torque  field.  The  charges! 
are  used  to  calculate  the  molecular  dipole  moment.  The  model  is  illustrated 
iin  Fig.  4.  Because  only  a  single  molecule  is  considered  here  (i.e.  no 
neighboring  molecules),  inertial  effects  associated  with  the  surrounding 
water  structure  is  very  crudely  represented  by  placing  hydrogen  atom 
(proton)  masses  at  those  point  charge  positions.  CMD  or  other  detailed 
theoretical  justification  does  not  exist  for  such  mass  placement. 

The  fixed  on  the  molecule  primed  coordinate  system  of  Fig.  4  is 
j  imbedded  in  the  laboratory  (x,  y,  z)  system  where  the  +z  direction  is  the 
{direction  of  shock  propagation  (see  Fig.  5).  The  calculation  proceeds 
by  using  Euler  angle  transformations  in  order  to  find  the  torque  about  the 
line  of  nodes.  The  torque  results  in  an  equation  of  motion 
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i where  the  sin  6  term  is  exact  (if  fn  is  exact)  and  the  second  term  is  a 
crude  restoring  force  with  m  being  a  proton  mass,  squ  the  oxygent-hydrogen 
intramolecular  distance  (1A),  k  is  the  hydrogen  bond  bending  force  constant] 
and  I  is  the  moment  of  inertia  of  the  model  molecule  about  the  line  of 
inodes  with  the  oxygen  atom  as  center. 

t 

The  shock  Induced  effective  force  (in  the  +z  direction)  on  each 
hydrogen  atom,  fn,  used  here  is  given  by 


fn  - 
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>  IT  is  the  peak  stress  amplitude,  Tg tHe^ffecVTrdnt* ^hietaess^ii^  isrthe  i 
pre-shocked  state  numberdenslty  of  '  water  moleciile‘sf  ~  and  a  Is  the  fraction- 
of  the  molecule  force  which  acts  on  a  hydrogen  atom  (a  «  1/18  for  H20), 

Eqs.  (2)  and  (3)  represent  a  three  dimensional  version  of  a  model  first 
attempted  by  Horie?* 

Equation  (2)  has  two  interesting  limiting  solutions.  One  solution  Is 
where  the  shock  induced  torque  is  j ust~balnWfe_d  by  therhydrogen  bond  bend¬ 
ing  restoring  force  so  that  5  ■  0.  The  second  solution  corresponds  to  the 
bond  broken  regime  and  thus  has  k  -  0. 

,  «« 

For  60  assumed  smaTl~~the  0*0  solution  can- be  shown  to  lead  to 


<cos0>  -  - 
Similarly  the  k 

i 

<COS0>  -  - 


Trlona§Hfn 

(L2)*ik  I  1! _ 

-  0  solution  can  be  put  In  the  form 

7Tto*at.i_ 

(48)**  I 


(4) 


(5) 


jin  arriving  at  Eqs.  (4)  and  (5)  the  angle  average  of  cos0  over  a  unit 
j  sphere  has  been  used.  Eqs.  (4)  and  (5)  give  the  average  shock  Induced 
| projection  of  the  modified  rigid  four-point-charge  along  the  direction  of 
i shock  propagation. 

i 

i 

Equations  (4)  and  (5)  are  employed  to  arrive  at  the  shock  Induced 
j  polarization  P 

P  -  npQ  <cos0>  (6) 

where  n  Is  the  number  density  of  molecules  which  contribute  to  the  observer 
shock  polarization  signals,  and  p©  is  the  four-point«-charge  molecule^ 
electric  dipole  igment  «  2.27  x  10“18  esu  cm  from  Fig.  4  with  |Vi|  - 
|v2|  -  1A  and  |V3|  -  |vA|  -  0.8  A).  Because  Ls  -  U8x8h  (with  the 

shock  front  rise  time),  Eqs.  (5)  and  (6)  predict  P  proportional  to  T8h 
for  t  -  Tsh. 

i  7  8 

|  In  turn  P  ls  used  to  calculate  the  prompt  pulse  current  denlsty  * 

!  j  .  p,  {Up  -  up)  (7) 

I  J  K  Lo 

where  K  Is  dielectric  constant,  lo  ls  the  sample  thickness,  and  up  the 
particle  velocity.  The  form  of  Eq.  (7)  is  appropriate  to  cgs  units. 
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where  e0  is  the  vacuum  permittivity  (rationalized  mks  units)  • 

POLARI Z  ATI  Otf^teTAr^NAbYSTS - - - — - 

@  Unbonded  (i.e.  k-0)  molecules. 


For  X r*  100  kbar  -cgs)— t- »-  t  4  ■  10  —  sec,  no  -  3.47  x  10^  cm~ 

Ls  -  Usxsh,  U8  «  5.05  x  105  em/see— (£ro»-Ref .  10),  and  I  -  3.50  x  10“*0 


gm  cm* 


>  -<cob&> 


-Ov8-2~y- 


-Cax 


which  is  equivalent  to  saying  that  all  of  the  contributing  shock  rotated 
molecules  are  almost  perfectly  aligned  with  the  direction  of  shock 
propagation  at  Tsh  =*  10*12  sec. 

-12 

For  V(0+)^»  0.2  Volt  one  has  agreement  between  the  xsh  *  10  sec 
voltage  rise  time  and  the  right  hand  side  of  Eq.  (8)  providing  that 


—  -  2.7  x  10' 


-3 


(10) 


The  Eq.  (10)  result  ls  within  a  factor  of  2  of  the  100  kbar  and  correspond^ 
ing  Hugonlot  temperature^  ionic  concentrations  obtained  by  using  the 
temperature  dependent  dissociation  constants  from  Eigen  and  MaeyerU  and 
the  empirical  conductivity  versus  hydrostatic  pressure  results  of  Holz- 
appel!2. 

The  Eq.  (10)  result  ls  also  within  a  factor  of  2  agreement  with  the 
0.5Z  ion  concentration  at  130  kbar  quoted  by  A,C.  Mitchell  el  al^  as 
a  consequence  of  behind  the  shock  front  electrical  conductivity  measure¬ 
ments.  It  is  also  within  a  factor  of  2  when  compared  with  the  ion  con¬ 
centrations  obtained  from  behind  the  shock  electrical  conductivity 
measurements  of  Hamann  and  Linton^  and  an  ion  mobility  of  3  x  10~3  cm2 
V-l  sec”*  (the  atmospheric  pressure  25°C  value)  from  Eigen  and  Maeyer*!. 

Applying  Lq  -  1mm,  and  (U8  -  up)  -  3.06  x  10^  cm/sec  to  Eq.  (7)  give^ 
0.2  Volt  for  the  voltage  predicted  across  a  50  G  load  by  the  experimental 
cell  of  0.5  cm^  cross  section  providing  that  K»Kq  and 
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“v  are  used.  Although  was  used  because  these  authors  do  not  know _ 

|  better.  It  should  be  pointed  out  that  in  so  far  as  It  produces  charges 
which  can  be  separated  In  an  electric  field,  a  high  degree  of  lonlclty 
|  will  Increase  the  dielectric  constant.  An  Increased  K  would  result  in 
j  better  agreement  between  Eqs.  (10)  and  (11). 

1 2 

^  *  Using  Tgh  »  10"  sec  and  the  above  k»0  formulae  yields  the  same  very 

good  agreement  between  theory  and  experiment  down  to  35  kbar.  Near  35  kbar 
the  slope  of  the  polarization  voltage  versus  pressure  changes  (see  Fig.  3)  ► 
►  i  The  implication  is  a  changing  shock  polarization  mechanism. 

}  (b)  Hydrogen  bonded  molecules* _ 

i 

j  Applying  the  previously  used  numerical  values  at  10  kbar,  and  a  k 

value!5  of  3,78  x  10-13  erg/radlan^,  to  Eq.  (5)  gives _ 


[  <cos0>  *■  -  1.5  x  10  .  (12) 

i 

I  In  turn  Eq.  (12)  yields,  via  Eq.  (7)  and  the  very  reasonable  bonded 

|  regime  assumption  that  on  the  average  each  shock  front  enveloped  water 
!  molecule  contributes  equally  to  the  shock  polarization  current  density 
i  (i.e.  n  *  nD) ,  for  the  predicted  prompt  peak  polarization  voltage  across 
[  a  50ft  load 

j 

S  V(50ft)  =  15  m  Volts.  (13) 

j 

)  K-Kq  has  been  assumed. 

I 

For  an  extrapolation  of  the  Fig*  3  experimental  data  to  6mV  at 
!  10  kbar,  the  left  hand  side  of  Eq.  (8)  becomes  6  x  10^  Volts/ sec  while 
the  right  hand  side  is  (from  Eq.  (12))  8.7  x  10^  Volts/sec.  It  is  truly 
amazing  that  the  very  crude  bending  angle  restoring  force  model  utilized 
in  Eq.  (2)  can  give  such  excellent  agreement  between  theory  and  experiment 

DISCUSSION 

Employing  the  well  known  shock  velocity  values^  andTsh=*  10~*^  sec 
one  finds  Ls  (100  kbar)  =  50A.  That  shock  front  thickness  value  is 
!  within  a  factor  of  2  agreement  with  the  computer  molecular  dynamics 
!  results* 6  of  Tsai  and  Trevino  for  shocks  in  dense  monatomic  liquids.  50A 
j  is  consistent  with  the  5.8  kbar  laser  reflectivity  value  upper  limit  of 
>  400A. 


Because  of  the  lack  of  experimental  data  in  the  low  pressure  region 
(i.e.  <20  kbars),  the  5.8  kbar  reflectivity  work  becomes  the  glue  which 
I  holds  the  work  over  the  entire  pressure  range  together.  It  can  be  shown 
'  t  hat  JthfeJfcfiOA jupperJUmi jL.pc.ex entn^Jve„uQk^^  _ taim. 
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*  being  the  explana t ion~f o  r-the-abservedlaw_pr  essure-p 

USSR  researchers*^"^  have  claimed  a  shock  front  related  viscosity, 
jU,  value  in  the  100  kbar  region  of  approximately  10^  p,  By  using 

0  ^  LspUs  (14) 

ifrom  Gilman^l,  which  should  be  valid  in  the  100  kbar  range  (p  is  mass 
[density),  one  calculates  n  ^  0*2  P  at  100  kbar  (from  Ls  -  50A). 

!  _ _ 

There  has  been  c  onsiderabLejils  cuss  ion  in -the  literature  relating  to 
the  USSR  value  of  10^P*  While  the  reasons  for  the  disagreement  are  not 
clear,  it  should  be  noted  that  Hamann  and  Linton^,  and  Hataann^,  have 
consistently  claimed  that  the  USSR  n  values  are  incorrect. 

- - ; - 

_ The JCshJ!LlJL_^aec^ndJAfc. corresponding-  claim  with  respect  to 

molecular  orientation  within  the  shock  front  (100  kbar  to  35  kbar  range) 
represents  a  major  advance  in  the  interpretation  of  shock  front  physics* 
Second,  the  detail  and  accuracy  of  the  experiment  and  theory  for  the  shock 
front  optical  reflectivity  serves  as  the  low  pressure  anchor  for  the 
entire  program  reported  on  here*  Further,  it  should  be  noted  that  the 
400A  upper  limit  reflectivity  result  limits  the  many  USSR  claims  that  such 
j reflectivity  measurements  can  predict  molecular  dimension  shock  front 
jthlcknesses;  the  thickness  is  too  small  in  comparison  to  the  optical 
(wavelengths  to  yield  sufficient  theoretical  accuracy  for  those  claims. 

|  That  the  bond  broken  regime  dominates  the  polarization  signals  at 
I  high  pressures  is  not  surprising;  Graham^  has  advanced  the  concept  that 
1  shock  Induced  bond  scission  is  the  agent  responsible  for  shock  polarization 
land  increased  electrical  conductivity  in  polymeric  solids*  At  the  moment, 
however,  it  does  not  appear  that  bond  breaking  is  the  polarization 
producing  agent  in  the  10  kbar  region* 
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|  Figure  3,  Averaged  experimental  data.  The  predicted  slope  and  amplitude 
i  at  10  kbar  is  for  net  alignment  of  hydrogen  bonded  dipoles • 
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